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ABSTRACT 

HST WFPC2 images are presented that span the inner ~19 kpc diameter of the edge-on galaxy NGC 
3079; they are combined with optical, emission-line imaging spectrophotometry and VLA images of 
radio polarization vectors and rotation measures. Ionized gas filaments within 9-kpc diameter project up 
to 3 kpc above the disk, with the brightest forming the 1 kpc diameter superbubble. They are often 
resolved into strands 0."3 (25 pc) wide, which emerge from the nuclear CO ring as five distinct streams 
with large velocities and velocity dispersions (FWHM «450 km s" 1 ). The brightest stream emits w 10% 
of the superbubble Ha flux and extends for 250 pc along the axis of the VLBI radio jet to one corner 
of the base of the superbubble. The other four streams are not connected to the jet, instead curving 
up to the vertical ss 0.6 kpc above the galaxy disk, then dispersing as a spray of droplets each with 
~ 10 3 \/JMq of ionized gas (the volume filling factor />3x 10~ 3 ). Shredded clumps of disk gas form a 
similar structure in hydrodynamical models of a galaxy-scale wind. The pattern of magnetic fields and the 
gaseous kinematics also suggest that a wind of mechanical luminosity L w w 10 43 ergs s" 1 has stagnated in 
the galaxy disk at a radius of ~ 800 pc, has flared to larger radii with increasing height as the balancing 
ISM pressure reduces above the disk, and has entrained dense clouds into a "mushroom vortex" above 
the disk. Ha emissivity of the filaments limits densities to n e > 4.3/~ 1,/2 cm" 3 , hence kinetic energy and 
momentum to (0.4-5) x 10 55 y/f ergs and (1.6-6) x 10 \pf dyne s, respectively; the ranges result from 
uncertain space velocities. A prominent star-forming complex elsewhere in the galaxy shows a striking 
spray of linear filaments that extend for hundreds of parsecs to end in unresolved "bullets." 

Subject headings: galaxies: individual (NGC 3079) — galaxies: active — galaxies: jets — galaxies: kinematics and 
dynamics — galaxies: ISM — galaxies: magnetic fields 
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1. INTRODUCTION 

An energetic outflow can inflate a "superbubble" to 
a size that exceeds the scale height of the interstellar 
medium (ISM). If the bubble ruptures, ejecta can al- 
ter chemical abundances across the galaxy. The su- 
perbubble is comprised of various shock fronts in the 
ISM and the outflow ejecta. The ISM and ejecta have 
several phases with different density, temperature, em- 
bedded magnetic fields, and dust distributions, and the 
shocks have different properties depending upon the 
phase through which they are propagating. The flow 
must therefore be observed across many wavebands 
from the radio to hard X-rays, but is often seen most 
clearly by imaging spectrophotometry of optical emis- 
sion lines. Spatial resolutions and sensitivities of cur- 
rent instruments suffice to diagnose gaseous conditions 
for only the nearest and brightest superbubbles (e.g. 
M82, Shopbell & Bland-Hawthorn 1998; NGC 253, 
Forbes et al. 2000; NGC 4438, Kenney & Yale 2000; 
Kenney et al. 1995) where the angular scale is large 
enough to constrain dynamical models. 

One of the clearest examples is the superbubble 
emerging from the disk of the edge-on SB-pec (de 
Vaucouleurs et al. 1991) galaxy NGC 3079; galaxy 
properties are summarized in Table 1 . The X-shaped 
radio structure on the nucleus discovered by de Bruyn 
(1977) was subsequently resolved into a figure-8 pat- 
tern by Duric & Seaquist (1988, hereafter DS) in the 
radio continuum, and by Ford et al. (1986) in optical 
line emission. The nucleus contains the most luminous 
known H2O megamaser (Henkel et al. 1984) which 
Trotter et al. (1998) resolved spatially, and has large 
extinction, so conditions there have been established 
from IR emission lines and radio-frequency spectra (Ir- 
win & Seaquist 1988; Irwin & Sofue 1992; Baan & 
Irwin 1995) which show outflow signatures. 

Heckman, Armus, & Miley (1990) and Filippenko 
& Sargent (1992, hereafter Paper I) obtained long- 
slit spectra of the ionized gas associated with the 
superbubble, and found large velocity gradients and 
emission-line widths in the patterns expected for an 
expanding structure. In Veilleux et al. (1994, hereafter 

2 SOAR Telescope Project, NOAO, Tucson AZ 85726-6732 
'Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. These obser- 
vations are associated with proposal ID GO-6674. 
3 Cottrell Scholar of the Research Corporation 



Paper II) we used Fabry-Perot spectrophotometry to 
map the ionized superbubble, showing that its pattern 
of radial velocities is consistent with an outflow whose 
space velocities increase from the nucleus. 

Because of its considerable extinction (Ptak 1997; 
Sofue et al 2000), the nuclear power source manifests 
itself primarily by reprocessed radiation. How the 
superbubble filaments are powered is thus uncertain. 
Their total ionized mass and kinetic energy (KE) are 
scaled by gas densities whose values are poorly con- 
strained due to blending in ground-based spectra (of 
the red [S II] doublet flux ratio) and images (to estab- 
lish the emitting volume). DS modeled the radio struc- 
ture as a wind-blown bubble from the active galactic 
nucleus (AGN), and Hawarden et al. (1995) concluded 
that an AGN wind is responsible largely because hot 
dust characteristic of a starburst is absent; Israel et al. 
(1998) subsequently found hot dust out to a radius of 
300 pc. The structure appears to be in a "breakout" 
stage, with the driving wind expanding freely into a 
halo of diffuse X-ray emission (Pietsch, Trinchieri, & 
Vogler 1998, hereafter PTV). An AGN in NGC 3079 
opens the possibility of a more collimated outflow, and 
indeed a jet has been proposed (Irwin & Seaquist 1988; 
Trotter et al. 1998; Sawada-Satoh et al. 2000) to ex- 
plain several aligned, VLBI-scale radio knots along 
with the maser kinematics (Sawada-Satoh et al. 2000). 
Irwin & Seaquist (1988) argued that a jet alone may 
inflate the superbubble. 

If the wind arises from an optically extinguished 
AGN, it provides a nearby laboratory for study of 
processes thought to have operated in dust-shrouded 
proto-quasars. If the wind is associated with a star- 
burst, the processes operating in IR ultra-luminous, 
distant galaxies can be studied in detail. To address 
these issues, in §2 we present Hubble Space Telescope 
(HST) WFPC2 and unpublished VLA 4 radio contin- 
uum images of the galaxy disk and superbubble. In §3 
we derive properties of the magnetic fields, and show 
that ionized filaments are often resolved transversely 
so that we can constrain their gas densities and fill- 
ing factor from recombination fluxes. Some filaments 
with large velocity dispersions align with the VLBI- 
scale axis, extending the directly detected influence of 
the jet to 250 pc radius. In §4 we discuss the dynamical 
state and power source of the superbubble. We outline 

4 The VLA is a facility of the National Radio Astronomy Observa- 
tory. The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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some future observations in §5, then summarize our 
results and conclusions in §6. 

2. OBSERVATIONS AND REDUCTIONS 
2.1. HST WFPC2 Images 

For sharpest images in exposures of reasonable du- 
ration, we spanned the inner parts of the galaxy with 
two chips of the WFC but "dithered" the exposures 
between two points separated by 2 "50 (confirmed to 
±0"02 by cross-correlating image pairs); see Table 
2. We imaged the galaxy on two rather than three 
WFC detectors to maximize coverage perpendicular to 
— rather than along — the plane of the galaxy disk. 
The PC chip covered the outskirts of NGC 3079, and 
did not obtain useful data. Our filters passed the I- 
band (F814W) and [N II]AA6548,6583+Ha (F658N) 
emission-line complex. The two I-band exposures 
were sufficiently short that we did not need to ob- 
tain multiple images to identify cosmic -ray hits. Four 
narrow-band F658N images were made, two at each 
dither position. 

"Warm pixels" were treated using the appropriate 
dark frame that became available several weeks af- 
ter our exposures. The dithered frame gave us infor- 
mation at these locations, so we interpolated linearly 
each frame to the position of its partner then patched 
in the affected pixels. The handful of residual cosmic 
rays and warm pixels in the combined dither images 
were removed manually. Because only a single contin- 
uum image was obtained at each dither position, cos- 
mic rays were removed by first interpolating linearly 
each frame to the position of its partner, then selecting 
the minimum pixel value at each point from this image 
and its uninterpolated partner. Six counts were added 
to the uninterpolated frame (an offset that was deter- 
mined empirically) so that the interpolated frame con- 
tributed only at the locations of cosmic rays or warm 
pixels in the uninterpolated frame. Less than 0.1% of 
the pixels in each frame were in this state, very few 
across the superbubble. 

After cleaning artifacts and scaling to a common 
total exposure, we interleaved the frames at the two 
dither positions onto a grid of 0"0707 pixels. In the 
continuum image 90% of a star's light is encircled 
within 0."38 diameter (0."226 ± 0"003 FWHM Mof- 
fat function), as expected for a two-point dither on the 
WFC chips at this wavelength. We scaled the contin- 
uum image fluxes to ensure a "line"-only image after 
subtraction. To place the I-band image on an absolute 



flux scale, we used the color term for an SBc galaxy 
given in Table 6.2 of the WFPC2 Instrument Hand- 
book. 5 

We obtained additional images from the HST 
archive (P.I. H. Falcke), see Table 2. Signal/noise ra- 
tios are smaller than those of our images because the 
galaxy nucleus was placed on the PC CCD and shorter 
exposures were used. However, the I-band archive im- 
age has slightly higher resolution on the nucleus than 
our dithered WFC images, and places the N half of 
galaxy on one of the WFC CCD's, complementing our 
coverage. Composite images are shown in Fig. 1. 

In the F658N-F814W image, the flux contributions 
of the three emission-lines [N II]AA6548,6583+Ha 
passed by the narrower filter vary across the disk and 
especially across the superbubble. In Veilleux, Ce- 
cil, & Bland-Hawthorn (1995, hereafter Paper III) and 
Veilleux et al. (1999, hereafter Paper IV) we presented 
flux maps and velocity fields, respectively, for ionized 
gas outside the superbubble. We showed in Paper II 
that gas in the interval 1 < R < 2.5 kpc has average line 
widths of w 150 km s" 1 , dropping to 100 km s" 1 far- 
ther out. We used these values to reconstruct emission- 
line profiles of the [N II]AA6548,6583+Ha complex at 
ground-based resolution (1" FWHM), and then multi- 
plied the result by the F658N filter profile to find the 
average flux fraction contributed by Ha at each point. 
Velocity gradients across the superbubble are large and 
vary in a complex fashion, see Fig. 2. In Paper I we 
had decomposed the profile of each spectral line of the 
[N II]AA6548,6583+Ha complex into 3 Gaussian sub- 
systems within rectangular regions 0"57 highx2."3 
wide. Using this decomposition, we synthesized the 
complexes, multiplied them by the F658N profile, and 
inserted the resulting Ha weights into our weighted 
image. Fig. 3 shows the distribution of weights. Cor- 
rections exceed those for the disk: at the base of the 
bubble only 25% of the line flux through filter F658N 
comes from Ha (the spike at left in Fig. 3). 

Next, we used the WFPC2 Exposure Calcula- 
tor 6 to find the Ha fluxes of H II regions whose 
[N II]A6583/Ha flux ratios were measured in our 
ground-based data. From the sum of counts at these 
regions, we placed Ha fluxes in our interleaved, line- 
only frame on an absolute scale. The major change 
over Fig. lb is to darken the lower half of the bubble 
where [N II] A6583/Ha > 1 (Paper I); see Fig. 4a. 

5 www.stsci.edu/instraments/wfpc2AVfpc2_hand/wfpc2_handbook.html 

6 www.stsci.edu/wfpc2/etc.html 
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2.2. VLA Flux, Polarization, and Rotation Mea- 
sure Images 

Radio emission traces the outflow on all observed 
scales, and many images have been made following 
the discovery of minor-axis radio lobes by Duric et 
al. (1983). Using VLA images at 1.4 and 4.9 GHz, 
they constrained the rotation measure to the E bubble, 
50 < R m < 300 rad m~ 2 . But polarization studies are 
tough at 1.4 GHz because the 30% polarization seen 
at 4.9 GHz has been reduced to <3% by Faraday rota- 
tion. From the many datasets of NGC 3079 in the VLA 
archive, we therefore extracted unpublished datasets 
made by Dr. S. Caganoff and collaborators which are 
centered at 8.1 and 4.4 GHz (Fig. 5 and Fig. 6a). De- 
tails of the observations are given in Table 2. 

Standard procedures under AIPS reduced the data 
at four frequencies. We then fit xi v ) = x(0) + A 2 /? m 
to each pixel to map the transverse electric field (E±) 
X(0) and the R m . Uniform weights produced a beam 
of 1.65"x 1.59" (HPBW) along PA = 13.3°. Figs. 5 
and 6 are insensitive to features >30". Fig. 6b shows 
that the 3.8 cm radio emission is < 3.4 x 10" 5 Jan- 
sky arcsec" 2 inside the superbubble. 

3. EMPIRICAL RESULTS 

3.1. Dust Plumes and the Pattern of Scattered 
Light in the Galaxy Bulge 

Much of the "peanut bulge" discussed in Paper 
IV may arise from scattered nuclear or inner-disk 
starlight. Israel et al. (1998) argue from their near- 
IR broad-band images that 20-30% of the bulge light 
must be scattered. The pattern we see in the I-band 
image (Fig. la) is bilobal and straddles the superbub- 
ble. The outflow appears to have depleted the cen- 
tral cylinder of dust, and indeed, several dust plumes 
arch out from the bottom half of the superbubble to 
fall back onto the galaxy disk. The increased dust 
load in the presumed toroidal lobe enhances scattered 
starlight. Scattering in the foreground halo also seems 
to "soften" the appearance of dust filaments near the 
nucleus and in the background disk to produce the 
striking impression of "depth" in e.g. Fig. 4a. 

We quantified dust extinction by comparing the I- 
band surface brightness measured toward a dust fea- 
ture Sdc.i with its local background St g j (e.g. Howk & 
Savage 1997). Foreground starlight reduces absorp- 
tion by Aj = -2.5\og(Sdcj/Sbgj)- Defining x as the 
fraction of starlight emitted in front of the feature, 



Sdcj/Sbgj = x + (l-x)e~ T ' where 77 approximates the 
maximum extinction optical depth through the dust. 
We ignored light forward scattered along our line of 
sight, and also assumed that dust has small depth com- 
pared to stars so it acts as an external screen. We aver- 
aged 5 -pixel wide intensity cuts through a prominent 
dust arc in Fig. 4a at the edge of the superbubble. 

A single color cannot constrain tightly both x and 
77. By bounding x (0.4 < x < 0.8), we found that 
0.27 < Aj < 1.6 mag. (0.48 < A v < 2.8 mag, assum- 
ing R = 3 for the ratio of absolute to relative extinc- 
tion). Multiplying N H = 1.7 x 10 21 A V cm" 2 by 1.37 
to convert H to total mass (He + typical metal frac- 
tion), and treating the dust arc as a cylinder of radius 
0."3 and length 3", we found that the total mass is 
(1-4.7) x 10 5 Mq , larger if dust was destroyed during 
cloud acceleration. Despite this high mass, the dusty 
feature shows no line emission, underscoring the need 
to map the outflow in many ISM phases. 

3.2. Dereddening 

The J-K color map of Israel et al. (1998) shows en- 
hanced reddening W of the nucleus. This is also ap- 
parent in Fig. lb as an over subtraction of the I-band 
image from the F658N exposure, which forms the dark 
zone at the base of the superbubble. Fig. 4a shows that 
the ionized filaments of the superbubble originate from 
an elliptical region s=a 7" x 3" in extent, with patchy ob- 
scuration. On this scale there is a ring between 1 ."25 
(Sofue & Irwin 1992, from CO) and 3. "5 (Meaburn et 
al. 1998, from molecular H and radio continuum) radii 
that is rotating at s=a 330 km s -1 (Israel et al. 1998). Is- 
rael et al. (1998) find H 2 (v = 1 -» 0)5(1) emission of 
w 10 7 Lq to a height of w 100 pc and out to 300 pc 
near the galaxy disk plane. 

In Paper II we used long-slit spectra of the Ha and 
H/3 lines (corrected for Balmer absorption) to deter- 
mine an average reddening optical depth of t(Ho) w 5 
at the nucleus, w 1 in the bottom third of the superbub- 
ble, and zero beyond; the value on the nucleus agrees 
with that found by Israel et al. (1998) across the inner 
6" x 2". Accordingly, to deredden the filaments, we 
bracketed the extinction gradients with height by us- 
ing two linear trends such that A v = 2.75 at the base 
of the bubble, declining to A v = 1.75 at heights be- 
tween 3 and 4", and to no reddening by 8" and 10.5" 
height, respectively. The superbubble is projected on 
diffuse, ionized emission which is a mixture of true 
vertical structure and the outskirts of the background 
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galaxy disk. As we summed the flux of each filament, 
we therefore subtracted the average vertical gradient 
outside the superbubble. 

Isolated in this fashion, the superbubble Ha flux 
sums to w 10~ 13 ergs s" 1 cm" 2 with compact filaments 
emitting slightly more than half, the jet (§3.10) 10%, 
and diffuse emission 1/3, of the total. We found an Ha 
luminosity of the filaments of (3.5-3.9) x 10 39 ergs s" 1 
(depending on the reddening), comparable to what we 
derived from the ground-based spectra (Papers I & II). 

3.3. Constraints on Ionized Densities and Masses 

Fig. 8 shows the Ha luminosity function of fila- 
ments labeled in Fig. 7, assuming the smaller redden- 
ing discussed above. We estimated a minimum volume 
of each filament by assuming a cylinder with depth 
equal to its thinnest dimension (measured at the con- 
tour of 50% peak emission) on the sky; many fila- 
ments are barely resolved in one dimension — «0"25 
FWHM thick after accounting for instrument blur — 
and most at the top of the superbubble are unresolved. 
Assuming Case-B recombination conditions at 10 4 K 
and n p w 0.9n e , the electron density is n e = 1.8 x 
10 12 y / LHa/Vf cm -3 for gaseous filling factor / and 
emitting volume V. The flux-averaged, ionized density 
of the superbubble filaments is n e = 5.7/- 1/2 cm" 3 ; Fig. 
8 shows the population distribution and the distribution 
of ionized mass which sums to 1 .4 x \Q h \/J Mq . We 
used Monte Carlo techniques to quantify the uncertain- 
ties of these derived parameters. We assumed that the 
Ha luminosity is known to ±25% (half of this uncer- 
tainty comes from the reddening ambiguity, half from 
the distance to NGC 3079) and the volume to ±50% 
(which encompasses the range of filament widths that 
are resolved in the Ha image). Table 3 lists mean and 
extreme values of the density and mass of each fila- 
ment under these assumptions, and whether the fila- 
ment is resolved. 

In Paper II we found no evidence from the red [S II] 
doublet flux ratios that the ionized gas exceeds the 
low-density limit n e > 100 cm" 3 across the superbub- 
ble. This and our new HST results now constrain the 
gaseous filling factor / to > 0.3% of the gas volume. A 
typical filament therefore has ionized mass > 300Mq ; 
filaments near the top of the superbubble are smaller 
and fainter, and we find have ionized mass > 60Mq . 
If filaments are fragments of molecular clouds lofted 
from the galaxy disk (see e.g. the models of Schiano 
1986; Suchkov et al. 1994), then this constrains only 



the mass of their outer ionized sheath. 

3.4. Superbubble Filaments 

Fig. 4 shows that four prominent bundles of ionized 
filaments emerge from the CO ring, and appear to be 
spaced «90° apart in azimuth. They are close together 
in (perhaps intertwined?) bundles until they reach a 
height of 550 pc, at which point each bundle separates 
and the strands, otherwise unchanged in appearance, 
make sharp angles on the sky while still trending up- 
ward. Fig. 6b shows that at greater heights strands are 
replaced by a spray of loops and arches that appear to 
"drip" toward the galaxy disk; the apparent connection 
between knots and arches may be spurious because 
of limited spatial resolution. Fig. 4b shows that the 
soft X-ray emission peaks along the drips. Beyond, 
isolated small clouds are evident > 1 .5 kpc above the 
galaxy disk. Beyond the clouds, radio emission contin- 
ues outward as a "cap" to 28" (2.3 kpc) radius. At still 
larger radii, soft X-ray emission fades into the ROSAT 
background (PTV). 

3.5. Constraints on the Magnetic Field 

Fig. 5c shows that the E radio lobe has a well- 
ordered magnetic field, indicated both by its high po- 
larization (20-30% at 6 and 3.8 cm wavelengths) and 
by the trend in rotation measure R m . While R m is 
nearly constant along the polar axis, perpendicular it 
changes sign from <~ +50 rad m~ 2 on the inside to 
<~ -50 rad " 2 toward the outer edge. The line-of-sight 
component of the magnetic field reverses dramat- 
ically whereas the transverse field B± (= E± + tt/2) 
aligns with the outflow, particularly along the N edge. 

The E lobe appears as a loop or limb-brightened 
shell; similar radio structures are seen in other inclined 
galaxies (e.g. NGC 2992, Wehrle & Morris 1987). 
Whether it forms part of a shell, closed or open loop 
or bubble, anchored or otherwise to the central disk is 
an important clue to its origin. Unfortunately, Fig. 5 
shows that we have information only for the sides of 
the radio lobe. The connection with the disk and the 
cap of the radio lobe (Fig. 5a) is too faint to derive po- 
larized intensities. However, the R m inversion - inside 
to outside - is compelling and, to our knowledge, has 
never before been observed. 

The form of R m excludes a field centered on the 
lobe because the projected field would integrate to 
zero. Our fits to R m follow a A 2 dependence closely, 
all but ruling out a Faraday-rotating medium within 
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the synchrotron source (Sokoloff et al. 1998). Be- 
cause polarizations are sometimes close to the max- 
imum possible value (75%, Burn 1966), little depo- 
larization can have occurred. Indeed, the slow spa- 
tial/spectral behaviour of R m compared to the beam 
rules out beamwidth/bandwidth effects. Thus, the 
Faraday screen lies predominantly in front of the syn- 
chrotron source and is partly resolved. The mean 
R m w 50 rad m" 2 is easily generated and constrains the 
density of the screen. Duric et al. (1983) derive a mean 
field of 20 fiG for the E lobe. We derive a lower limit 
onB| by equating \R m \ w 50 rad m" 2 with the limit on 
Ha emission measure in the lobe (E m < 1 cm" 6 pc). If 
the depth of the Faraday screen is L ~ 10 pc, because 
R m = 0.81B||« e L, we derive Bn > 15 [iG. The required 
gas density is n e < 0.3 cm" 3 . 

At the wind speed postulated in §4.2.1, it is unlikely 
that magnetic energy density dominates the thermo- 
kinetic energy of the flow. The wind may simply 
carry along the magnetic field (Klein, Wielebinski, & 
Morsi 1988). Unless the compressed field is much 
higher than implied by minimum energy arguments, 
the above numbers are fairly restrictive and suggest 
that the radio lobe beyond the superbubble should emit 
detectable Ha emission at a flux of 10~ 18 ergs cm" 2 s" 1 
arcsec" 2 , w 20% of our current 2cr threshold (Paper II). 

In light of the optical data, we suspect that the po- 
larized emission arises from a limb-brightened incom- 
plete shell (partial bubble) that is tilted slightly to the 
line of sight. The radio bubble is expanding into a 
magnetized, hot (PTV), low-density halo. The radio 
continuum is presumably limb-brightened where the 
field is highly compressed. In §4.2.3 we provide a 
qualitative model to explain the R m inversion. 

3.6. Optical Emission from the Counterbubble 

In Fig. 1 1 of Paper II we traced line emission from 
high-velocity gas between projected radii of 10" and 
15" (panels c & d) W of the nucleus. Here Fig. lb 
shows line-emitting regions with [N II]A653/Ha flux 
ratios like those in the E superbubble. Spectra of IR 
recombination-line flux ratios would constrain the red- 
dening, hence energetics of the ionized counterbubble. 

3.7. Morphology of Disk H II Regions 

Properties of the star-forming, large-scale galaxy 
disk will be considered in a later paper. However, ener- 
getic, non-planar motions are of immediate relevance. 

A notable feature of the H II regions in NGC 3079 



is their tendency to "blister open" at their tops, see es- 
pecially the outer ring in the right-hand half of Fig. lb. 
H I arc "C" of Irwin & Seaquist (1990) lies above this 
region. As discussed in §4.4 of Paper III, such struc- 
tures suggest vertical motion of gas up "galactic chim- 
neys" which also allow photons to escape from the disk 
to energize the prominent diffuse ionized medium. 

A striking active region is marked near the left-hand 
side of Fig. lb, where Fig. la shows a high density 
of resolved stars and dust filaments. This complex 
lies well beyond the maximum radius of the stellar 
bar, is a local flux peak in the 245 GHz continuum 
(8 mJy over 11" diameter), lies below H I arc "D" of 
Irwin & Seaquist (1990) which requires formation en- 
ergies (0.5-2) x 10 54 ergs, and has X-ray luminosity 
w 5.4 x 10 38 ergs s" 1 in the ROSAT HRI band (source 
H14 of PTV). Fig. 9 shows that over a diameter of 15", 
roughly a dozen linear filaments protrude for 0.6 kpc. 
They have broad line widths (w200 km s" 1 corrected 
for instrumental resolution), but mean velocities equal 
to that of the rotating galaxy (Paper IV for the ion- 
ized gas, and Irwin & Seaquist (1991) for the H I) at 
their projected radius. Fig. 9 also shows unresolved 
knots of line emission at the ends of many of these fila- 
ments. The three brightest knots have Ha luminosities 
of (1.5 ±0.1) x 10 37 ergs s" 1 assuming no reddening. 
Two features labeled "?" in the figure have velocity 
FWHM w200 km s" 1 in our Fabry-Perot spectra (Fig. 
4c of Paper III), and sharp edges on the side directed 
away from the source of the linear filaments. These 
structures resemble superficially Galactic Herbig-Haro 
objects, but are far more energetic. 

3.8. Velocity Field of Ionized Filaments in the Su- 
perbubble 

3.8.1. Previous Results 

As mentioned in §3.2, emission line profiles in 
Fig. 2 have a component from gas in the background 
galaxy that is separable from energetic gas in the su- 
perbubble by its modest deviation from the galaxy 
systemic velocity 7 and by the small values of its 
[N II]A6583/Ha flux ratios (which are typical of H 
II regions, see Fig. 3 of Paper II). The spectral pro- 
file from a 0"57 highx2."3 wide box often spans sev- 

7 As discussed in §3.2.1 of Paper IV, the systemic velocity of NGC 
3079 depends on the waveband of observation. Values 1150±25 
km s -1 result from the combined effects of a disk bar and possi- 
ble warp, and differences in optical depth and orientation to various 
galaxy components. We use v sys = 1 150 km s" 1 . 
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eral hundred km s , arising from "turbulent" motions 
within a filament or gradients that are unresolved at 
ground resolution. An example is at the bottom left- 
hand corner of the superbubble, where filaments 45, 
49, & 61 have opposite velocities (±150 km s" 1 ) de- 
spite being adjacent on the sky. 

In Paper II we showed that the superbubble veloc- 
ity field is a partial "Doppler ellipsoid," with the cen- 
troids of emission-line components spanning radial ve- 
locities -1050 to +550 km s -1 relative to v sys , and hav- 
ing maximum velocity splitting 1250 km s" 1 near the 
mid-axis of the superbubble. Velocities do not return 
to small values at the top of the superbubble as they 
do in spatially resolved spectra of planetary nebulae; 
instead, line emission fades away as velocities reduce 
modestly from extreme values. In Paper II we showed 
that this pattern is quite well reproduced by a nuclear 
outflow wherein filament space velocities accelerate to 
3700 km s" 1 at r = 1 kpc as v{r) oc r a , with a ~2.5, 
while motion vectors swing toward the sky plane as 
the bubble cap is approached. The main deficiency is 
in the N half of the superbubble, where the model pre- 
dicts redshifts from v., yj of 550 km s" 1 vs. 200 km s" 1 
observed for the presumed backside gas. 

3.8.2. Refinement Motivated by Present Data 

The HST line image shows that the filaments are 
organized into four well-defined vertical bundles, each 
fortuitously spanning the 1" resolution of our ground- 
based spectra. Although we cannot assign an unam- 
biguous value to a filament, we can bound velocities 
in each bundle with the radial velocities mapped in 
Fig. 4 of Paper II. We find that most filaments have 
blueshifted motions relative to v sys , with largest radial 
velocities in the spray of arched filaments near the top 
of the bubble. Anticipating the dynamics to be dis- 
cussed in §4.2, it is plausible that these are "mush- 
rooming" out of the bubble with a substantial motion 
along our sightline. So we refined the field used in Pa- 
per II, now assuming expansion mostly parallel to the 
disk plane at set height in the bubble (Fig. 10), with 
space velocity comparable to the observed velocity of 
the blue wing along the mid-axis of the bubble. At left 
in Fig. 7 we show the adopted space velocity at each 
height. The expansion velocity that we infer at the 
base of the bubble is consistent with the values derived 
there by Hawarden et al. (1995) from H2 transitions 
and mapped by Israel et al. (1998) (see also Fig. 11). 
Together, this model and the nucleus-centered pattern 
used in Paper II should bracket the true motions. 



3.9. Bulk Energies and Momenta of the Ionized 
Filaments 

Table 3 lists maximum and minimum space veloc- 
ities of each filament as derived from the projections 
of Paper II and §3.8.2, respectively. Total KE and 
momentum of ionized gas in the superbubble sum to 
(2.9-0.26) x 10 55 V/ergs and (3.4-0.74) x 10 47 V7 
dyne s, respectively; values for each filament are also 
shown in Table 3. We retain our estimate from Paper 
II of "turbulent" energy — 3.1 x \0 54 /n e ergs for the 
blue-wing component. This is now an upper limit, be- 
cause we now resolve some of the line broadening into 
discrete filaments each with a bulk KE. 

3.10. Jet Emission Near the Nucleus 

Fig. 1 1 shows that the brightest stream of ionized 
filaments projects near the S boundary of the super- 
bubble base, and coincides with the axis of the pu- 
tative VLBI-scale jet mapped at 8 GHz by Trotter et 
al. (1998) and Sawada-Satoh et al. (2000). In the 
magnified inserts in Fig. 1 this feature has Ha sur- 
face brightness «5x 10" 15 ergs s" 1 cm" 2 arcsec -2 co- 
incident with a narrow plume of 1=16 mag arcsec -2 
continuum light that points to the nucleus. Using the 
HST line-only image Fig. lb to define a mask of the 
jet emission, we summed spectra from our ground- 
based Fabry-Perot datacube to map velocities (see Fig. 
1 lb), then averaged spectra at two intervals along the 
jet (Fig. 11c); the mean velocity of the two summed 
spectra is redshifted by 125 km s" 1 relative to v sys . 
Elsewhere the corrected FWHM is 200 km s" 1 and 
[N II]/Ha ss 1 .2. In the two summed spectra we sub- 
tracted this component to isolate a broader jet-specific 
feature whose Ha flux we then deblended from the ad- 
jacent [N II]AA6548,6563 doublet. This component 
has centroid blueshifted 125 km s -1 relative to v sys . Fit 
results are reported in Table 4; the KE is a lower limit 
assuming that the broader feature arises entirely from 
expansion around the jet that projects all motion onto 
our line of sight. Actually, gas kinematics constrains 
jet orientation and space velocity poorly because gas 
is ionized at radiative shocks within the thermally un- 
stable cocoon/ISM interface. 

4. DISCUSSION 
4.1. Nuclear Jet 

Irwin & Seaquist (1988) suggested that the super- 
bubble could have been blown by a precessing VLBI- 
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scale jet; Baan & Irwin (1995) elaborated on this sce- 
nario. Irwin & Sofue (1992) interpreted extensions 
of CO(J=l-0)-emitting gas into the base of the super- 
bubble as resulting from a jet/ISM cloud interaction. 
Using pressure confinement arguments and the mini- 
mum momentum flow, they estimated that the putative 
jet has velocity > 10 4 km s" 1 and power > 4 x 10 40 
ergs s" 1 , respectively. Over the 1 Myr dynamical age 
of the superbubble (Paper II) the KE exceeds 10 54 ergs, 
which suffices to power the CO outflow. 

Irwin & Seaquist (1988) resolved the nuclear emis- 
sion into a flat-spectrum radio core with parsec-scale 
collimation along P.A. 120°. Trotter et al. (1998) fur- 
ther resolved this structure into four clumps with radio 
spectra inverted between 5 and 8 GHz, and steep be- 
tween 8 and 22 GHz; in this sense the nuclear sources 
resemble a low-power analog of powerful radio galax- 
ies. Sawada-Satoh et al. (2000) found that the spectral 
indices and fluxes of several of these knots changed 
considerably in 22 months and that knot motions were 
consistent with expansion from the nucleus at an ap- 
parent velocity of <~ 0.16c. The four clumps align pre- 
cisely (see §3.10 and Fig. 1 1) with the brightest stream 
of filaments at the base of the superbubble. Three re- 
sults suggest that the elongated, ionized feature is en- 
ergetic (not ambient gas photoionized by the AGN): 

1. Fig. 4b shows that X-rays peak off-nucleus. 
ROSAT was pointed to « ±2" accuracy on this 
galaxy (PTV), about the displacement from the 
nucleus in the Figure. Ptak (1997) modelled the 
ASCA+ROSAT PSPC spectrum to show that the 
nucleus is quite obscured ((1-2.3) x 10 21 cm" 2 ), 
so the jet may be generating these X-rays. 

2. Velocities of both ionized (§3.10) and molecular 
gas (Israel et al. 1998) peak in this extranuclear 
region. Velocities of order V S 4oo = 400 km s" 1 
will produce shock emission L(Ha) s hock = 1 .5 x 
10 34 noAV 3 400 ergs s" 1 (A is the shock surface 
area in pc 2 in an ambient ISM of density no 
cm" 3 ). The luminosity of filament #3 in Table 3 
arises from pure shock ionization if no ~ 4 cm" 3 . 

3. Here optical emission-line profiles are broad (<~ 
440 km s" 1 FWHM, §4.1), as we would expect 
from gas shock-heated as a jet impacts the inner 
boundary of the molecular cavity in the disk. In- 
deed, Israel et al. (1998) see bright H2 emission 
and hot dust T w 900 ± 100 K nearby. 



Rossi et al. (2000) showed by numerical simula- 
tion that a jet propagating through the ISM in a spi- 
ral galaxy radiates on average 0.1-0.5% of its KE 
as ionizing photons at each interaction site, with ion- 
izing luminosity \5L H p over a wide range of jet ve- 
locities and ISM densities. From Table 3 the ionizing 
luminosity 1.7 x 10 39 ergs s" 1 can therefore inject the 
KE w 5 x 10 52 /0.005 = 10 55 ergs of the superbubble 
over its dynamical age if the jet has operated with high 
enough duty cycle. 

We can constrain the age of the radio bubble inde- 
pendent of the optical outflow. If we assume the bubble 
was inflated by the jet but then lifted by buoyancy, the 
buoyancy timescale is T b <~ 2RyJ r/GM{R) where r is 
the cavity radius and R is the distance to the center of 
the galaxy. This assumes that the surrounding gas has 
to fall around the cavity at speed v <~ y / GM(R)r/R 2 
and that the cavity must move <~ 2r before gas needs 
to be replaced. Our picture is of a light bubble which 
is temporarily overpressured relative to its surround- 
ings and thus rises and expands as the ambient density 
dropped. The timescale for buoyant rise always ex- 
ceeds that for reaching local pressure equilibrium. 

We take the cavity as 1 kpc diameter at 5 kpc from 
nucleus. We derive M(R) from the Galaxy mass model, 
M(R) ~ 10 10 fl kpc (V max /220) 2 M Q such that the mass 
goes as R (in kpc) beyond 3 kpc (Freeman 1996) as- 
suming a spherical potential. The last term involves 
Vmax which is the deprojected rotation velocity on the 
fiat part of the rotation, curve. In Paper IV we 
find for NGC 3079 that V max = 250 km s" 1 such that 
M(5 kpc) - 6 x 10 10 M Q , so T h ~ 1 x 10 7 yrs. 

To form the whole superbubble, the jet would have 
had to precess from its current orientation at the S cor- 
ner. However, the four vertical bundles of filaments 
that constitute the superbubble and which all break at 
the same height above the disk are not explained by 
jet precession or bouyancy. To inflate this structure we 
now examine a wide-angle wind as proposed by DS. 

4.2. Dynamics of the Wind-Blown Superbubble 

As an outflow inflates a superbubble, the various 
gas phases organize into a characteristic pattern of 
shock and ionization fronts, as described by e.g. Schi- 
ano (1985, 1986) and simulated by e.g. Strickland & 
Stevens (2000): from largest to smallest radii in a ma- 
ture superbubble there is undisturbed ISM, perhaps a 
photoionized precursor, a dense shell of compressed 
ISM that forms a standing "ring" shock, the con- 
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tact discontinuity between the two fluids, hot shocked 
wind, a wind shock, and free-blowing wind near the 
sites of energy injection. Several such structures may 
be evident, depending on the density structure of the 
galaxy ISM and on the detailed history of the power 
source. The hot, shocked wind cools by lateral ex- 
pansion and by mixing with the ambient ISM at the 
Rayleigh-Taylor (RT)- and Kelvin-Helmholtz-unstable 
interfaces between the two fluids. 

If the wind is sufficiently energetic or long-lived to 
"breakout", then its flow becomes a strong function of 
polar angle seen from the nucleus: free flow up to a 
critical angle wherein its ram pressure becomes com- 
parable to the thermal pressure in the diffuse compo- 
nent of the ISM; at larger angles a standing bow shock 
in the galaxy disk decelerates and deflects wind around 
a region of undisturbed ISM gas. The wind subse- 
quently reaccelerates into the halo. In the following 
subsections, we will identify many of these structures 
in our data. Table 5 lists dynamical parameters derived 
in Paper II, scalings we will use in §4.3. 

4.2.1. ISM "Ring" Shock 

Except for the volume occupied by the superbubble 
filaments, the concave bowl of a pre-blowout super- 
bubble must be filled largely with shocked wind (zone 
B in the nomenclature of Weaver et al. 1977). The free- 
wind zone A is on a much smaller scale, presumably 
within the cavity of the molecular ring in NGC 3079. 
If the wind has not encountered enough background 
gas to halt its free expansion, it blows out and then the 
shocked wind is found only near the outer (ISM) shock 
in a ring near the galaxy disk. 

We can infer from Fig. lib that the wind must 
have a high enough luminosity to have pushed the 
ring (ISM) shock to w 8" = 900 pc radius near the 
disk plane. The models of Schiano (1985) show 
that this occurs when L w > 10 42 ergs s" 1 for reason- 
able ISM properties for NGC 3079, and the result- 
ing flow can reach the galaxy halo if it persists. The 
smallest wind luminosity for blowout occurs when 
the diffuse ISM (i.e., neglecting dense clouds) swept 
up by the wind becomes comparable to the wind 
mass, at time f* by Schiano (1985). For an expo- 
nential variation of ambient gas density e~ R l a , we 
have v w t* = a when L w = 6.2 x lO^a^n-iv 3 , ergs s" 1 , 
with a the scale height in kpc, n_i the gas density 
in 0.1 cm" 3 , and vg the wind velocity in 10 9 cm s" 1 . 
The equal-pressure stagnation/stand-off radius in the 



galaxy plane Rep = aft = O^^L^/vg^P-io kpc, where 
ft 2 = L w /(2-KPisMV w a 2 ), a combination of wind and 
ISM parameters. At blowout, Schiano (1985) has ar- 
gued that the contact discontinuity between shocked 
wind and ambient ISM lies near the ISM ring shock, 
a conjecture supported by simulations (e.g. Strickland 
& Stevens 2000). Pressure balance then sets the angle 
ip between the radial wind direction R = \J r 2 + z 2 and 
the surface normal to the contact discontinuity 

x 2 ~t~ y 2 

cos 2 -0(x,y) = —^-h(x,y), 

where x = r/a,y = z/a, and h(x,y) is the pressure distri- 
bution. The shape of the contact discontinuity follows 
by solving the quadratic equation 

(y) 2 [x 2 - (x2+ ^ (x ^V 2x,y + [/- ( ^ 2+ ^^ ) ]=o 

for y' = dy/dx and then numerically integrating an il- 
parameterized model from the initial value in the disk- 
plane y(x = Rep) = 0. 

We considered several forms for h(x,y): exponen- 
tials in both x and y, exponential only in x (y constant 
over the height of the superbubble), and power-law 
only in x (again, y constant). Fig. 12 scales these forms 
to match the observed shape of the outer shock. The 
H I and diffuse ionized gas scale-heights of NGC 3079 
are uncertain but may be w0.5 kpc (Irwin & Seaquist 
1991, §3.1 of paper II, respectively) to 1 kpc, where- 
upon the filaments extend to y/a = 4-2. Evidently the 
calculation is not very constrained, but the simulations 
of Suchkov et al. (1994) to be discussed shortly do 
match our images with an exponential disk and spher- 
ical densities for the stellar spheroid and halo. 

4.2.2. Comparison to Hydrodynamical Models 

The HST images strengthen our assertion in Paper 

II that the superbubble has "blown out" of the galaxy 
disk: the linear X-filaments discussed in Papers II and 

III are now shown to join the galaxy disk at radii of 
?y 8", forming the "concave bowl" that we attribute 
to the standoff outer disk shock/contact-discontinuity. 
Other features also suggest a disk shock: 

1. Schiano (1985) estimated the evolving shape of 
the wind-driven complex, see his Fig. 8. For 
wind parameters logL w = 43,logv w = 9, and 
\ogPjsM = -10 (cgs units and values similar to 
those inferred by DS except for being twice their 
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wind velocity), the inner (wind) shock lies at 
two scale height radii where we see the filamen- 
tary shell, while the outer (ISM) ring shock sits 
at « 3 scale heights in the disk where we see the 
edges of the concave bowl. 

2. The structure has a large [N II]/Ha flux ratio that 
is correlated with line widths of 100-250 km s" 1 
(Paper III), consistent with shock emission. 

3. For central gas density 100 cm -3 and wind speed 
v w = 10 8 7 cm s" 1 , the wind model of DS has a 
central cavity at 0.5 kpc radius which expands 
at 300 km s" 1 . The radius is about half what we 
observe up the bubble (but similar to the bright 
green region at its base in Fig. 4), while the ex- 
pansion rate they expect is comparable to what 
in §3.8.2 we attributed to that of the back wall of 
the bowl. That gas reaches 550 km s" 1 relative 
to v sys . Hawarden et al. (1995) rescale the pa- 
rameters of the DS model to the higher densities 
expected in the molecular disk at a radius of 50 
pc, and obtain reasonable agreement with shock 
velocities inferred from their H2 spectra. 

Suchkov et al. (1994) and Strickland & Stevens 
(2000) simulated the appearances at X-ray and opti- 
cal temperatures of a superbubble at various stages of 
its wind-driven evolution. They parameterized models 
by temperatures, scale heights, masses, and densities 
of the galaxy disk and isothermal halo. Suchkov et al. 
did not include a stellar disk because the bulge gravity 
dominates; their model injects energy within a radius 
of 150 pc, comparable to the molecular disk in NGC 
3079 (Fig. 11). They discussed the time evolutions of 
the temperature and density distributions. 

To constrain parameters of NGC 3079 we used the 
X-ray results from PTV, who inferred from the spa- 
tial variation of ROSAT PSPC spectra (and the HRI 
image) of NGC 3079 a hot halo with density n e <~ 
7 x 10" 4 /" 1 / 2 cm" 3 ; depending on the average value 
of / in the ionized filaments, this n e is consistent with 
either of the Suchkov et al. "B"-models. In contrast, 
the "A" models have a large halo mass, and lead to 
flows that differ dramatically from what is observed in 
NGC 3079, viz. the relatively large mass above the 
disk causes the wind to breakout laterally across the 
disk, not what we observe. Mapping the putative hot 
halo awaits the increased sensitivity of the Newton X- 
ray Observatory (XMM); these data may provide evi- 
dence that the halo results from extensive star forma- 



tion in the galaxy disk that predated the superbubble. 

However, with current X-ray data we can say that 
the HST emission-line image more closely resembles 
the predictions of Suchkov et al. model B 1 rather than 
B2. Bl has lower disk gas density than B2 (no = 2 
vs. 12 cm" 3 ), is hotter (T disk = 6 x 10 4 vs. 1 x 10 4 K), 
and has ten times the halo density hence a higher disk 
vertical scale-height (160 vs. 60 pc). The hotter gas 
disk in B 1 seems more appropriate for a star-forming 
system like NGC 3079. The halo density of B2 is too 
tenuous to confine the wind. Bl produces a conical 
"nozzle" with length/diameter ratio comparable to that 
of the concave bowl in NGC 3079. In contrast, the 
tenuous halo of B2 is much less collimating. 

Peak velocities derived for disk-gas filaments en- 
trained in the wind are ~4300 km s" 1 for Bl vs. w 700 
km s -1 for B2 (Suchkov et al. 1994); these straddle ob- 
served values in the superbubble (Paper II). Fig. 13c 
shows velocities of Bl during breakout at 3.2 Myr 
(older than the superbubble in NGC 3079). At this 
stage of the model outflow, much of the space velocity 
near the top of the bubble is projected along our line 
of sight. The HST images do not suggest vortices until 
radii of w 10". The models predict that the vortices 
will eventually propagate down the walls and grow to 
disrupt the superbubble completely. 

Another feature of model Bl resembles the HST 
structures: because the halo is so tenuous, entrained 
disk material tends to shatter into individual fila- 
ments. Each cloud impedes the wind flow separately 
via a stand-off bowshock upstream in the wind. The 
shocked wind acts to confine hence preserve the cloud 
during its acceleration. Such behaviour persists in 
three-dimensional numerical simulations (e.g. Stone 
& Norman 1992). Overall, the B models reproduce 
many observed characteristics of this outflow. 

4.2.3. Wind Compression of the Magnetic Field 

The radio data also support a blowout vortex. The 
dramatic discontinuity in R m , high polarization lev- 
els, and aligned field may arise from the strong com- 
pressional shock of the wind. Models (Brandenburg, 
Moss, & Shukurov 1995, hereafter BMS) show that 
outflows are expected to entrain disk material and 
frozen-in magnetic field. While initially highly turbu- 
lent, the field becomes stretched and compressed along 
the flow. Here we associate the positive R m regions 
with the updraft material, see Fig. 14b. 

BMS show that as the entrained material cools, it 
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rains back toward the disk outside the flow thereby re- 
versing the field. This behavior is seen in the hydrody- 
namical simulations discussed in §4.2.2. By associat- 
ing sinking gas with regions of negative R m , we obtain 
a plausible way to reverse R m in the E lobe (Fig. 6c). 

BMS build their model on the observation that most 
halo fields appear to lie parallel to the Galactic disk 
(e.g. Beck et al. 1996). In our interpretation of Fig. 3 
in BMS, the halo field is not seen. When these mod- 
els are density weighted, only the central updraft and 
reverse flow are evident. The line of neutral magnetic 
field between regions of positive and negative R m is 
probably unresolved by our beam, and need not depo- 
larize significantly. 

An interesting variant on the BMS scheme is to 
associate the negative R m material with compressed, 
magnetized halo material (Fig. 14a). Global magnetic 
fields in galaxies are often well ordered and axisym- 
metric (Beck et al. 1996), and can extend many kilo- 
parsecs into the halo. Dynamo models are summarized 
by Zel'dovich, Ruzmaikin, & Sokoloff (1983, here- 
after ZRS). Disk and halo fields can have odd (dipole) 
or even (quadrupole) parity, and may even be of dif- 
ferent polarity. This forms the basis for our alternative 
model for the field reversal in Fig. 6c. Again the wind 
drags its own field, but now runs into a quadrupole par- 
ity halo field with polarity opposite to the wind (e.g. 
Brandenburg et al. 1993). The sign of the rotation mea- 
sure over the W lobe may establish if we are probing 
the compressed halo field or entrained material in out- 
flows. Either outcome would be of interest. 

A more exotic possibility (Yusaf-Zadeh, & Morris 
1987) is that the radio loop arises from a twisted field 
that protrudes from the disk, much like a solar promi- 
nence. This topology could explain the R m inversion. 
But, unlike a prominence, the loop cannot produce the 
observed synchrotron by winding up the field at the 
base of the loop. Magnetic energy in the twisted field 
increases at the rate W ~ B 2 v ( f ) A/&Tr, where A is the 
footprint area and is the winding velocity. If we as- 
sociate 1 % of the loop's radio energy with the winding, 
for a 300 pc diameter footprint the required twist rate 
is more than a hundred times larger than the vorticity 
or shear (<~ 10 km s" 1 ) associated with rotating disks. 
The disk field is entrained and dragged with the highly 
energetic flow. The wind is expected to evolve faster 
than the dynamo (10 6-7 yr vs. 10 8-9 yr, respectively). 
The mean magnetic field in the wind seems high. But 
circumnuclear activity may provide the turbulent dy- 
namo that feeds the azimuthal field through differen- 



tial rotation. This can be converted to vertical field if 
there are azimuthal gradients in the disk (Sokoloff, & 
Shukurov 1990). 

4.2.4. Disruption of Superbubble Filaments 

Our images resemble the models of Schiano (1985), 
suggesting that the ionized filaments in the superbub- 
ble lie near the wind shock. In detail, this shock may 
take the form of separate bow shocks upstream of each 
dense filament. Bathed in hot shocked gas that flows 
by subsonically, the cool filaments will erode from 
thermo- and hydro-dynamical instabilities. However, 
filaments have survived long enough for the thermal- 
ized wind to expand for several kpc to form the stand- 
off shocked boundary of the concave bowl. What does 
this interpretation of the data say about the filaments? 

Table 3 gives densities of the ionized outer sheath of 
the filaments where column length N c i ou d > A 'sheath > 
3 x 10 19 cm" 2 , whereas the estimate from dust in §3.1 
applies only to non-line emitting structures. Because 
filaments have lofted to ml kpc height, we can esti- 
mate their minimum column to survive saturated ther- 
mal conduction (Cowie & McKee 1977). To resist 
evaporation (Krolik, McKee, & Tarter 1981) requires 
that N cloud > 5 x 10 21 cm- 2 (/?^ c /v 60 o) 6/7 7; 4 1/7 P e 6/ 9 7 , 
where P Cj 9 is the thermal cloud pressure in nT = 10 7 
cgs units, Rkpc is the height of the surviving filaments 
in kpc, and V6oo is their space velocity in units of 600 
km s" 1 . The embedded magnetic fields (§4.2.3) are too 
weak to alter these numbers much (Mineshige, Shi- 
bata, & Shapiro 1993). 

Suchkov et a l. (1994) show that RT instabilities 
grow as tRj = \J ^zp/2T:vl irc Ap, where p is the com- 
bined density of the ISM gas, and Ap w 0.01. These 
authors note that the wavelengths of perturbations 
which eventually disrupt the outer shell of a superbub- 
ble are of order the shell thickness. In NGC 3079, we 
see that the N shell wall is unresolved, < 0."2 FWHM 
so A < 17 pc. From the model velocity field of the 
bar-forced disk in Paper IV, we find a circular velocity 
of 100 km s" 1 at the radius of the superbubble. With 
these numbers, tgj ss 0.85 Myr. The dynamical time 
to inflate the bubble is somewhat greater (Table 5), so 
it is plausible that RT instabilities shear filaments into 
arches at the top of the superbubble. Because the ISM 
filaments are organized into discrete streams, they ap- 
pear never to have formed a dense shell; thus, the wind 
can indeed break out into the halo with ease. 
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4.3. Constraints on the Power Source 

Assuming the minimum value of / <~ 3 x 10~ 3 from 
§3.3, the power source must be able to inject KE 
> 1.4 x 10 53 ergs into the ionized superbubble; this is 
comparable to the KE in the starburst-driven outflow 
of M82 (Shopbell & Bland-Hawthorn 1998). Because 
the W radio lobe probably has a similar complex of op- 
tically obscured filaments, the KE of the ionized out- 
flow is comparable to that driven by the AGN in NGC 
1068 (w 4 x 10 53 ergs, Cecil, Bland, & Tully 1990). 

We showed in Paper II that the starburst age and de- 
pletion time for gas to form stars exceed considerably 
the age of the superbubble. In support of numerous 
supernova remnants within a nuclear starburst, Pedlar 
et al. (1996) showed that the steep/inverted-spectrum 
radio sources are embedded at 20 cm in nonthermal 
emission that extends to 4" along PA. 170°. Meaburn 
et al. (1998) showed that this structure coincides with 
the warped molecular disk found by Sofue & Irwin 
(1992) to have total cool gas 10 10 M o and w 1500 M Q 
of H 2 (v = 1 — > 0)5(1) gas at its center. Braine et al. 
(1997) found that dust emission at 1.2 mm within a di- 
ameter ~10" is much less centrally concentrated than 
either the CO or radio continuum emission; they used 
this emission to derive a low N(H2) / /co(2->i) conver- 
sion ratio, and hence a more modest mass of cold 
gas near the center, «2x 10 8 Mq (adjusted for 17.3 
Mpc distance). The center appears not to contain cool 
dust (T < 20 K), thereby resembling a low-luminosity 
version of an IR ultraluminous galaxy. Braine et al. 
(1997) therefore argued that stars are forming effi- 
ciently in this region, whereas Sofue et al (2000) find 
an ultra-high-density molecular core and maintain that 
its large differential rotation inhibits star formation. 

To match L w the required injection rate of KE from 
a nuclear starburst would be w 3 x 10 41 ergs s" 1 (Table 
5). We had noted in Paper II that the observed mid-IR 
flux is consistent with a starburst of this power, ass- 
suming no contribution from AGN emission that may 
also have been reprocessed by dust. Using the injec- 
tion rate that we derived from our new data in §3.9 and 
assuming an adiabatic superbubble, we therefore ob- 
tained 3 x 10 41 > 2.6 x 10 54 V7/^>«- Assuming that 
the outflow is no larger than the maximum observed 
velocity of 1050 km s" 1 and that our reduced estimate 
of the KE in the superbubble is correct within a factor 
of 4, this condition is satisfied even for / = 1 . 

If radiative losses are important, the bubble may 
still conserve momentum. Then the observed injec- 



tion rate > 7.4 x 10 46 y7 /td yn dyne can be compared to 
that predicted from the starburst model of Elson, Fall, 
& Freeman (1989) p* = 6 x 10 33 L,>,n dyne w 2 x 10 33 
dyne. Again, a starburst suffices (in some cases just 
barely) for all values of /. 

Our data constrain only the ionized outflow. En- 
trained neutral and molecular material (Israel et al. 
1998; Irwin & Sofue 1992) and the wind reduce the 
acceptable value of / for the ionized gas. If the higher 
space velocities of our Paper II model are favored over 
those in this paper, then a starburst wind will work only 
for very small values of /. 

5. FUTURE PROSPECTS 

5.1. Spectra of the Warm, Post-Shock Gas 

Its large angular size is well resolved by current X- 
ray telescopes, so the superbubble in NGC 3079 of- 
fers excellent prospects to study confined hot gas in 
a galaxy-wind/ISM interaction. Both PTV and Ptak 
(1997) have modeled the PSPC spectrum obtained by 
ROSAT; Ptak also included ASCA data (no spatial res- 
olution but higher spectral resolution). However, large 
contamination by the relatively brighter galaxy disk 
and halo made their models ambiguous. The Chandra 
X-ray Observatory (CXO) now greatly improves spa- 
tial resolution. Unfortunately, Strickland & Stevens 
(2000) have demonstrated by simulation that the su- 
perbubble dynamics are likely to remain poorly con- 
strained because gas that emits soft X-rays has small 
filling factor compared to the shocked wind (which at 
> 10 8 ' 5 K is too hot to see with present instruments). 

However, CXO can distinguish between shocks and 
cooling gas at various T near the jet. As discussed in 
§4.2.2, hot gas likely clumps to the w 1" resolution 
of CXO as it cools behind strong high-velocity (>400 
km s" 1 ) shocks to 10 4 K. Re gions with pre-shock den- 
sity no (in units of cm" 3 ) would then be separated by up 
to the cooling distance of 88V s 4 00 /«o pc = 1 x V S 5 00 /«o 
arcsec, where V S 5oo is the shock velocity in units of 
500 km s" 1 . The distance will be somewhat smaller 
because most of the shock is oblique, not normal. 

5.2. Internal Properties of the Ionized Filaments 

To better understand their dynamics, we must con- 
strain the ionized column lengths, hence total masses, 
of the filaments. The dust calculation in §3. 1 and ther- 
mal evaporation timescale estimated in §4.2.4 are use- 
ful first steps, but column lengths are more directly 
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constrained with maps of emission-line fluxes. The ve- 
locity field of ionized gas associated with the jet must 
also be constrained. Integral-field spectrometers on 8- 
m-class telescopes can begin this work. Images of the 
filaments with higher spatial resolution may soon be 
possible with a laser-assisted, adaptive optics system 
on these telescopes (NGC 3079 does not have a natu- 
ral guidestar). These images would further constrain 
/ and, together with constraints on the column length, 
can exclude decisively a starburst-driven wind on en- 
ergetic grounds. 

Mid-IR spectroscopy with SIRTF will be able to 
constrain the cutoffs of the IMF, and the spatial and 
temporal evolution of putative starburst activity in the 
nuclear region, as has been done with the ESA Infrared 
Space Observatory for M82 by Forster-Schreiber 
(2000). With photo- and shock-ionization models such 
spectra could establish if filaments are pressure con- 
fined, and whether they have sufficient dust load for 
efficient wind acceleration. With filament properties 
understood, the superbubble of NGC 3079 will be- 
come a laboratory to study the galaxy-wind driven 
processes thought to prevail at high-redshifts. 

6. SUMMARY 

We have analyzed VLA spectra and HST emission- 
line and I-band images of NGC 3079. Using the ob- 
served gas distribution and recombination fluxes at 
high spatial resolution, we have tightened constraints 
on the densities of ionized filaments that delineate the 
superbubble. The HST images and pattern of mag- 
netic fields determined from radio polarization and ro- 
tation measures suggest independently that gas at the 
upper boundary is lofting out of the bubble in a vortex. 
Numerical simulations show that a substantial fraction 
of the space velocities of these filaments would then 
project along our sightline. We therefore reassess the 
velocity field across the bubble and reduce substan- 
tially previous estimates (including ours in Paper II) of 
the minimum outflow KE and momentum. We find that 
a bright, kinematically disturbed filament at the base of 
the superbubble and a narrow, straight plume of contin- 
uum light align with the putative VLBI-scale radio jet, 
extending the signature of the collimated component 
of this outflow to a projected radius of ^250 pc. The 
rest of the superbubble is composed of four gaseous 
columns whose orientation perpendicular to the disk 
plane and strikingly similar appearances are not easily 
explained by precessing jets. Either starburst or AGN 



power sources remain viable on energetic grounds be- 
cause of uncertainties in the flux of reprocessed in- 
frared photons from the AGN. But either source will 
have more modest steady-state requirements based on 
our data and kinematical model. 
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Fig. 1 .— HST WFPC2 images of NGC 3079, showing the log of intensity, with P. A. -10° at top and ticks every 30" = 
2.5 kpc. The top parts are noisier because they come from shorter, archival exposures. The I-band image (a) has been 
unsharp-masked to increase the contrast and prominence of the bilobal pattern of scattered light and dust filaments in 
the galaxy bulge. The image of [N II]A6583+Ha line emission (b) shows narrow ionized filaments that rise above the 
galaxy disk to envelope the prominent superbubble. The region spanned by Fig. 9 is delineated. The nuclear region 
is enlarged in the inset panels, with the optically obscured nucleus marked with a "+"; the I-band panel was obtained 
with the PC CCD of WFPC2, the line image is from dithered WFC CCD frames. Note the plume marked in both 
panels, which aligns with the VLBI-scale radio jet. (c) Composite of the [N II]A6583+Ha emission-line and I-band 
images in red and green, respectively, and the VLA 3.8 cm continuum in blue. 
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Fig. 2. — Fabry-Perot spectra of the superbubble (from Paper II), summed along the indicated "long slits", and oriented 
with PA. 80° at top; ticks are every 4" vertically and every 500 km s" 1 horizontally. Complex velocities across the 
superbubble blend the lines in the [N II],Ha complex. The transmission profile of WFPC2 filter F658N is shown at 
bottom right, with dashed lines intersecting at 50% peak transmission. Arrowheads mark the galaxy systemic velocity 
for the two bright spectral lines. Note broad spectral lines at the base of the superbubble, labeled "jet" in panels (e) 
and (f). 
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Fractional Contribution of Ha to F658N Counts 

Fig. 3. — How Ha contributes to the total flux measured through WFPC2 filter F658N. 
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Fig. 4. — (a) A concave bowl is shown enveloping the superbubble in a region 50" wide; its base appears to be devoid 
of H II regions. I-band starlight is green, Ha red, and [N II] blue. Note that some dust plumes align with ionized 
filaments in the bottom third of the superbubble while others at left and right arc up and away from the bowl. The 
brightest plume at the base of the superbubble aligns with the VLBI-scale jet (Trotter et al. 1998) (see also the inserts 
in Fig. 1). (b) Contours of the ROSAT HRJ emission (values given in Fig. 4 of PTV) are shown atop the HST line 
image, with registration uncertainties of ±2". The box delineates the region in Fig. 11a, while the white dashed ellipse 
shows the orientation of the galaxy disk. 
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Fig. 5. — (a) Polarized 6 cm flux is shown as an inverted greyscale. (b) Same as (a) now for polarized 3.8 cm flux, (c) 
Derived rotation measures, greyscale spans ±100 radians m" 2 . (d) The box delineates the same region as in the other 
panels. Contours are at levels (90, 50, 30, 20, 5, 2.5, 1, 0.6, 0.55, 0.5, 0.45, 0.4, 0.35, 0.3, 0.25, 0.2, 0.15, 0.1, 0.05, 
-0.05) percent of the peak surface brightness of 96.2 mJy beam -1 at 3.8 cm. Lines show the magnetic field vectors. 
Ticks in (a)-(c) are every 2", N is at top, and each panel is 20" on a side. 
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Fig. 6. — (a) The 6 cm VLA image with ROSAT HRI contours across the central part. The radio beam is 1 "59 x 1 ."65. 
The region shown is 4' wide, and is oriented like Fig. 1. Same contours are shown as in Fig. 4b. (b) Line emission 
(green) is superimposed on the radio continuum (yellow) image in (a). Representative areas are labeled in units of 
3.7 x 1CT 4 Jy arcsec" 2 . Magnetic field vectors, obtained by fitting P.A. vs. A 2 , are drawn on the radio structure. The 
polarization of the bright knot at top is uncertain due to 2ir ambiguities. 
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Fig. 7. — Filaments in the HST line image are numbered for reference in the list of their physical properties, Table 
3. The orientation is that of Fig. 2, with P.A. 80° at top and ticks every 0."5. The vertical lines show the effective 
"slit width" and sampling of our imaging Fabry-Perot spectrophotometry (when binned as in Fig. 2). Within each 
"slit" the dashed line shows the variation of radial velocity in the range (right hand side) to -1000 (left) km s -1 
relative to galaxy systemic. These trends are used to assign first an average radial velocity to each filament (see text), 
then — assuming cylindrical symmetry — the minimum space velocity shown at the left-hand side of this figure. 
For comparison, the velocities used in Paper II are shown in parentheses. Filaments with negative (approaching) 
radial velocities are labeled in black or white, receding velocities are shown in red. We assume that extinction varies 
vertically across the superbubble, and show at bottom right the multiplicative factor used to deredden the Ha fluxes at 
each height for the range of reddenings consistent with the observed Balmer emission-line decrements (Paper II). 
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Fig. 8. — Distributions of the derived properties of the ionized filaments in the superbubble. Ha luminosity is in units 
of 10 38 ergs s" 1 , and assumes the smaller reddening correction (see text). Electron density is in units of f~ x l 2 cm" 3 . 
Ionized mass is in units of 10 4 ^/JMq , and assumes case-B recombination conditions at 10 4 K. 




Fig. 9. — Linear filaments marked with lines emanate from a major star-forming complex in this line-emission image. 
The area shown is the 40" wide boxed region near the left-hand edge of Fig. lb, and lies beyond the maximum radius 
of the stellar bar. The line at bottom spans 15"! .25 kpc. Two larger, apparently edge-brightened nebulae that may be 
associated with this outflow are marked with "?". A structure resembling a star-forming ionization front is circled. 



21 




Fig. 10. — Inferred velocity fields associated with the nuclear outflow. The orange dashed lines point to highly 
inclined ionized filaments that may form the stagnation boundaries of the nuclear wind as it blows through the rarefied 
component of the galaxy ISM. The green lines point to the dense component of the disk ISM which has been elevated 
to form the superbubble. The bubble has ruptured into vortices at its top (schematic motions are shown in yellow). 
The red line at its base shows the filament stream that aligns with the VLBI-scale jet (see Fig. 11). Circled at bottom 
in dashed green are filaments whose emission-line spectra show high velocities and non-H-II region excitation (Paper 
II) perhaps associated with the optically extinguished radio counterbubble. 
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11" = 925 pc 
Ionized filaments 






Fig. 11. — (a) The base of the superbubble in Ha emission near the nucleus, shown with the same orientation as Fig. 
2. The image was smoothed with a Gaussian kernel of 0"18 FWHM. The brightest filament stream trends toward the 
S (left) boundary of the superbubble, and coincides in projection with the putative VLBI-scale jet, which is shown 
magnified 250 times using the image of Trotter et al. (1998) (not all radio-jet knots are shown). The bottom of the 
image is obscured by a dusty molecular ring, (b) Shown is the space-velocity diagram of the "long-slit" spectrum 
extracted from our Fabry-Perot datacube along the jet in a 1 "7-wide box. The two most prominent spectral lines are 
labeled, while the white lines at lower left show the regions over which the jet emission-line profiles were extracted 
for display in panel (c). (c) Emission-line profiles averaged from the regions shown in (b). Each spectral bin is 0.57 A 
wide («33 km s" 1 ). Arrowheads indicate v S ys for the line species, with the wavelength scale increasing to the right. 
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Fig. 12. — The shape of the contact discontinuity /ISM shock for various density distributions of the ambient ISM is 
shown to scale with the superbubble filaments. The curves delineate the equal-pressure stagnation distance at different 
hieghts above the disk. The dimensionless parameter £1 = 0.5 for all curves, reasonable for this outflow (see text). The 
solid line plots the shape (see §4.2.1) for exponential density functions in both coordinates. The ISM density depends 
only on radius for the other lines: dotted has form x a with a = -3/2, dashed has form e ax with a = -1 .9. 




Fig. 13. — Densities, temperatures, and velocities from model Bl of Suchkov et al. (1994). The first two panels are 
shown at time 2 Myr, the third at 3.2 Myr (i.e. after gas blowout) for a model system with more energetic parameters 
than NGC 3079. For these simulations, image size is 2 x 4 kpc 2 , the range of log« is [-3.9, 0.8] in cm" 3 , the range of 
log T is [4,8] in K, and the outflow velocity reaches 3 100 km s" 1 . Note that much of the space velocity can be projected 
along our line of sight (the eyeball depicts our perspective on the superbubble). 
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(a) halofiefd + windfiefd 



\b) wind field with reversaF 




Fig. 14. — Two models which attempt to explain the inversion of rotation measures by invoking a field reversal in 
a wind that is inclined to our line of sight. In (a), regions of R m > arise from the inner surface of the compressed 
radio-emitting shell. The magnetized wind (B w ) expands into a diffuse magnetized halo (B q ) with quadrupole parity 
(and opposite polarity). Here, regions of R m < arise where the wind has compressed the halo medium. In (b), 
the dragged field reverses where the entrained medium cools out of the flow. The outflow produces R m > and the 
returning medium reverses the field to give R m < 0. 
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Table 1 

NGC 3079 — Basic information 



Parameter 


Region a 

(") 


Flux 


Luminosity 


Reference 




total 


10.45 mag 


3.1 x 10 10 


1 


X-ray 


core 




3.0 x 10 6 


2 


X-ray 


lobes 




- 8 x 10 6 


2 


X-ray 


3 disk sources 


each 


1.6 x 10 5 


2 


12 /tin 


total 


2.81 Jy 




3 


25 /zm 


total 


3.54 Jy 




3 


60 /im 


total 


52.8 Jy 




3 


100 urn 


total 


96.5.8 Jy 




3 


40-122 /im 


total 




4.9 x 10 10 


3 


350 fim 


-80 


10.0±1.9 Jy 




4 


450 fim 


-80 


2.1±0.9Jy 




4 


800 Aim 


53 


0.92±0.1 Jy 




5 


1220 ztm 


total 


0.5 Jy 




6 


25 /im 


11 


1.5 Jy 




6 


60 /im 


11 


22.3 Jy 




6 


100 fim 


11 


40.8 Jy 




6 


800 /im 


11 


0.35 Jy 




6 


1220 /tm 


11 


0.03 Jy 




6 


CO(l-O) 


21 


231 Kkms" 1 




6 


CO(2-l) 


11 


306 K km s" 1 




6 


HI 


total 




7.1xlO y M Q 


7 


distance 


17.3 ± 1.5 Mpc 






8 


pos. angle 


169±4° 






8 


inclination 


82±4° 






8 



a "total" = value for entire galaxy, the units of most other entries in the column 
are diameters in arc-seconds. 



References.— (1) de Vaucouleurs et al. (1991); (2) PTV; (3) Young et al. 
(1989); (4) Eales, Wynn- Williams, & Duncan (1989); (5) Hawarden et al. (1995); 
(6) Braine et al. (1997); (7) Irwin & Seaquist (1991), with h = 0.75; (8) Paper IV 
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Table 2 
Images of NGC 3079 





WFPC2 




Very Large Array 


Time on Source at a 


Filter 


Source b 


Exposures c 


Configuration 


Date 


4.5 GHz 


8.06 GHz 






(s) 






(s) 


(s) 


F814W 


11/98 


400, 400 


B 


1/13/92 


6660 


6915 




Archive 


70, 70 


C 


12/9/90 


13110 


14670 


F658N 


11/98 


1400+2700 












11/98 


1700+2700 












Archive 


400, 400, 400 










F547M d 


Archive 


160, 160 











a Two 50 Mhz bands were used around each central frequency. The flux standards were 
3C 48 and 3C 286 (Dreher, Carilli, & Perley 1987). 

b Our HST exposures were made in November 1998; archive exposures were made in 
March 1999 and were reprocessed with on-the-fly recalibrations in November 2000. 

c "+" indicates the exposures at two dithering points. 

d This image has low signal/noise ratio on the superbubble, so was not scrutinized. As 
of November 2000 the public HST archive also contained a STIS UV exposure of 1320 s 
at central wavelength A 147 nm (P.I. D. Calzetti). It spans the base of the superbubble and 
nucleus, but shows only a few bright H II regions W of the nucleus. 
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Table 3 

Physical Properties of the Superbubble Ionized Filaments 



ID F Ha X 10 15 


Space Velocity" L Ha X 1(T 37 


Volume 


b ne c 


MassxlO" 3 d 


KExl0~ 51 e 


Momentum X 10 -45 e 


(erg s _1 crrr 2 ) 


(kms -1 ) (ergs -1 ) 


(PC -3 ) 


(l/vT'cm- 3 ) 


(V/M ) 


(Vferg) 


(v 7 / dyne s) 



1 


0.15 


75,500 


0.54 


<2999 


14.0 [10-22] 


1.2 [0.7-1.7] 


0.07, 3.13 


0.02, 


0.12 


2 


2.59 


200,500 


9.46 


<4()223 


15.9 [12-25] 


18.8 [12.3-25.5] 


7.46, 46.60 


0.75, 


1.87 


3 f 


9.6 


125,500 


32 


372146 


9.6 [5-14] 


150 [100-300] 


11.6, 1850 


1.8, 


7.4 


4 


0.25 


-200,-500 


0.915 


<5565 


13.3 [10-20] 


2.2 [1.4-2.0] 


0.87, 5.45 


0.09, 


0.22 


6 


4.56 


125,500 


16.83 


82576 


14.9 [11-23] 


35.8 [22.0-48.4] 


5.56, 89.02 


0.89, 


3.56 


7 


2.07 


-200,-500 


7.56 


<26102 


17.7 [14-27] 


13.5 [8.8-18.5] 


5.36, 33.51 


0.54, 


1.34 


9 


1.72 


-200,-500 


6.28 


<22677 


17.3 [14-27] 


11.4 [7.0-15.8] 


4.55, 28.45 


0.46, 


1.14 


11 


15.1 


-350,-500 


56.9 


381657 


12.5 [9-19] 


139.4 [88.0-193.7] 


169.88,346.69 


9.71, 


13.87 


12 


1.18 


-200,-500 


4.45 


<43214 


10.3 [8-16] 


13.1 [8.8-18.5] 


5.22, 32.63 


0.52, 


1.30 


13 


2.97 


-250,-750 


11.23 


44074 


16.4 [14-26] 


21.0 [15.0-29.1] 


13.08,117.74 


1.05, 


3.14 


15 


3.58 


175,750 


13.8 


94132 


12.3 [9-19] 


33.8 [22.0-48.4] 


10.30,189.18 


1.18, 


5.06 


16 


1.08 


-350,-750 


4.18 


<28241 


12.3 [10-19] 


10.2 [7.0-14.1] 


12.44,57.12 


0.71, 


1.53 


19 


1.03 


-350,-750 


3.96 


<19251 


14.5 [11-23] 


8.2 [5.3-11.4] 


9.97, 45.80 


0.57, 


1.23 


20 


0.98 


-200,-750 


3.77 


34232 


10.7 [8-17] 


10.7 [7.0-15.0] 


4.23, 59.55 


0.42, 


1.58 


23 


1.00 


0,750 


3.87 


<30806 


11.4 [9-18] 


10.2 [7.0-14.1] 


0.00, 0.00 


0.00, 


0.00 


24 


2.03 


-250,-750 


7.81 


<44501 


13.4 [10-20] 


17.4 [12.3-24.6] 


10.84, 97.53 


0.87, 


2.61 


25 


0.27 


-350,-1000 


1.05 


<8556 


11.2 [9-18] 


2.8 [1.9-4.0] 


3.43, 28.03 


0.19, 


0.55 


26 


1.23 


-250,-1000 


5.10 


99270 


7.0 [6-11] 


20.3 [13.2-29.9] 


12.64,202.27 


1.01, 


4.05 


28 


1.60 


275,1000 


6.63 


<57761 


10.5 [8-17] 


17.7 [13.2-24.6] 


13.31,176.01 


0.97, 


3.52 


33 


1.18 


-25,-1000 


4.89 


<59474 


8.9 [7-15] 


15.4 [10.6-22.0] 


0.10,154.94 


0.08, 


3.17 


34 


1.69 


-100,-1000 


6.995 


79585 


9.2 [7-15] 


21.3 [15.0-29.9] 


2.12,212.16 


0.42, 


4.23 


36 


0.39 


-250,-1000 


1.60 


< 16686 


9.5 [8-16] 


4.7 [3.5-7.0] 


2.90, 46.34 


0.23, 


0.92 


37 


2.38 


-250,-1000 


9.80 


<162161 


7.6 [6-12] 


36.1 [26.4-52.8] 


22.43,358.90 


1.80, 


7.18 


38 


0.47 


-50,-1000 


1.96 


<34658 


7.4 [6-12] 


7.5 [5.3-10.6] 


0.18,73.95 


0.08, 


1.58 


40 


0.15 


300,1000 


0.63 


<6843 


9.3 [7-15] 


1.8 [1.3-2.6] 


1.66, 18.39 


0.11, 


0.38 


41 


0.44 


-250,-1000 


1.81 


62899 


5.2 [4-9] 


9.7 [6.2-14.1] 


6.02, 96.34 


0.48, 


1.94 


42 


1.23 


0,1000 


5.09 


62039 


8.9 [7-15] 


16.1 [11.4-22.9] 


0.00, 0.00 


0.00, 


0.00 


43 


0.38 


375,1000 


1.55 


< 16259 


9.5 [7-16] 


4.6 [3.1-6.2] 


6.40, 45.51 


0.34, 


0.92 


44 


0.08 


400,1000 


0.34 


<3426 


9.8 [8-16] 


1.0 [0.7-14.1] 


1.54,9.63 


0.08, 


0.20 


45 


0.87 


-150,-1000 


3.59 


90280 


6.1 [5-10] 


16.3 [10.6-23.8] 


3.64,161.98 


0.48, 


3.23 


46 


0.59 


100,1000 


2.45 


<29520 


9.0 [7-15] 


7.7 [5.3-11.4] 


0.77, 77.47 


0.16, 


1.58 


47 


0.34 


-250,-1250 


1.395 


41935 


5.7 [5-9] 


7.0 [4.8-10.6] 


4.32,108.06 


0.34, 


1.72 


48 


0.22 




0.89 


<27600 


5.6 [4-9] 


4.5 [3-6.5] 








49 


1.90 


-125,-1250 


7.84 


128789 


7.6 [6-12] 


28.8 [20.2-40.5] 


4.47,447.21 


0.71, 


7.13 


50 


1.13 


-100,-1250 


4.656 


< 109105 


6.4 [5-10] 


20.4 [14.1-29.1] 


2.03,317.75 


0.40, 


5.06 


51 


2.54 


500,1250 


10.47 


151040 


8.2 [7-14] 


36.0 [25.5-49.3] 


89.52,559.51 


3.58, 


8.96 


52 


0.70 


-375,-1250 


2.89 


<32945 


9.2 [7-15] 


8.8 [5.3-13.2] 


12.32,136.84 


0.66, 


2.20 


53 


0.63 




2.60 


< 375970 


2.6 [2-4.5] 


28.3 [20-40] 








55 


0.30 




1.25 


<34230 


5.9 [5-10] 


5.9 [4-8.5] 








58 


0.09 


250,1250 


0.377 


< 19685 


4.3 [3-7] 


2.5 [1.8-3.5] 


1.53, 38.29 


0.12, 


0.62 


60 


0.27 


-250,-1250 


1.11 


< 17546 


7.7 [6-12] 


4.0 [2.6-5.7] 


2.46,61.62 


0.19, 


0.97 


61 


1.38 


150,1250 


5.76 


127077 


6.6 [5-11] 


24.4 [16.7-35.2] 


5.46,379.04 


0.73, 


6.09 


66 


0.20 


-375,-1250 


0.796 


<50491 


4.0 [3-6] 


5.8 [3.5-8.4] 


8.13,90.28 


0.43, 


1.44 


67 


0.42 


-400,-1250 


1.685 


<42362 


6.2 [5-10] 


7.7 [5.3-11.4] 


12.32,120.36 


0.62, 


1.93 


70 


1.04 


-550,-1250 


4.19 


105253 


6.2 [5-10] 


19.3 [13.2-27.3] 


58.00,299.57 


2.11, 


4.80 


71 


0.08 


-375,-1250 


0.307 


<22250 


3.6 [3-6] 


2.4 [1.3-3.5] 


3.33, 36.97 


0.18, 


0.59 


73 


0.20 


125,1250 


0.822 


31233 


5.1 [4-8] 


4.7 [2.6-7.0] 


0.72, 72.19 


0.11, 


1.14 


75 


0.06 




0.233 


<26500 


2.9 [2.5-5] 


2.3 [1.6-3.3] 








76 


0.13 




0.519 


< 13250 


6.2 [5-10] 


2.4 [1.6-3.3] 








77 


2.07 


-100,-1250 


8.16 


200671 


6.4 [5-10] 


37.2 [26.4-52.8] 


3.70,577.72 


0.74, 


9.24 


78 


0.05 




0.20 


<9390 


4.6 [3.5-7.5] 


1.3 [0.8-1.8] 








79 


0.14 




0.553 


<35330 


3.9 [3-6.5] 


4.1 [2.5-6] 








80 


0.19 


-400,-1500 


0.773 


<43214 


4.2 [3-7] 


5.3 [3.5-7.9] 


8.41,118.23 


0.42, 


1.58 
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TABLE 3 — Continued 



ID F Ha x 10 15 


Space Velocity 11 


L Ha x 1CT 37 


Volume 


b n e c 


MassxlO" 3 d 


KExl0~ 51 e 


Momentum X 1CT 45 c 


(erg s~' cm 4 ) 


(kms- 1 ) 


(erg s -1 ) 


(pc- 3 ) 


(1/V7cm- 3 ) 


(v?m ) 


(Werg) 


(VJ dyne s) 



82 


0.78 


200,1500 


2.768 


< 47066 


8.1 [6-12] 


11.1 [7.0-15.0] 


4.41,248.09 


0.44, 


3.30 


83 


0.78 


-125,-1500 


2.778 


103548 


5.5 [4-9] 


16.6 [10.6-22.0] 


2.57,370.17 


0.41, 


4.97 


84 


1.54 


200,1500 


5.504 


91993 


8.2 [6-12] 


21.9 [13.2-29.1] 


8.72,490.73 


0.87, 


6.54 


89 


0.09 




0.313 


<20430 


4.1 [3-6.5] 


2.5 [1.7-3.5] 








90 


0.11 




0.403 


<21530 


4.6 [3-7] 


2.9 [2-4] 








92 


0.85 


-650,-1500 


3.02 


< 69742 


6.9 [5-10] 


14.2 [8.8-19.4] 


59.56,317.16 


1.83, 


4.23 


93 


0.06 


150,1500 


0.227 


<9416 


5.2 [3-8] 


1.4 [0.9-1.9] 


0.32, 31.69 


0.04, 


0.44 


94 


0.65 


550,1750 


2.306 


44501 


7.6 [6-11] 


9.9 [6.2-13.2] 


29.66,300.26 


1.08, 


3.45 


96 


0.12 




0.432 


<2982() 


4.0 [3-6] 


3.5 [2.4-4.5] 








97 


0.14 




0.513 


< 32020 


4.2 [3-6.5] 


3.9 [2.5-5.5] 








98 


0.14 




0.499 


<33130 


4.1 [3-6] 


4.0 [2.5-5.5] 








99 


0.04 


-425,-2000 


0.154 


< 12834 


3.6 [3-6] 


1.4 [0.9-1.8] 


2.53, 55.95 


0.12, 


0.58 


101 


0.37 


-550,-2000 


1.327 


<28667 


7.2 [6-11] 


6.0 [3.5-7.9] 


18.01,238.17 


0.65, 


2.37 


103 


1.06 


225,2000 


3.779 


92845 


6.7 [5-10] 


18.2 [10.6-24.6] 


9.17,724.79 


0.82, 


7.28 


104 


0.12 


250,2000 


0.412 


<29954 


3.9 [3-6] 


3.4 [2.2-4.4] 


2.13,136.35 


0.17, 


1.34 


105 


0.05 




0.183 


< 18770 


3.3 [2.5-5] 


1.8 [1.2-2.3] 








106 


0.12 


-625,-2000 


0.438 


<11121 


6.6 [5-10] 


2.1 [1.3-2.9] 


8.20, 84.02 


0.26, 


0.85 


107 


0.04 




0.154 


< 13800 


3.5 [2.5-5.5] 


1.4 [0.9-1.9] 








109 


0.18 


-500,-2000 


0.631 


< 19685 


6.0 [5-9] 


3.4 [2.2-4.4] 


8.54,136.63 


0.34, 


1.37 


110 


0.10 


-900,-2000 


0.370 


<1 1121 


6.0 [5-9] 


1.9 [1.3-2.6] 


15.60, 77.04 


0.34, 


0.76 


112 


0.31 


-625,-2000 


1.09 


39789 


5.6 [4-9] 


6.4 [3.5-8.8] 


24.97,255.66 


0.80, 


2.56 


114 


0.03 




0.117 


<17120 


2.7 [2-4.5] 


1.4 [0.9-1.8] 








116 


0.001 




0.004 


<7730 


0.8 [0.5-1.2] 


0.2 [0.1-0.25] 








117 


1.36 


-900,-2000 


4.86 


117234 


6.8 [5-10] 


23.2 [15.0-29.1] 


187.21,924.51 


4.16, 


9.25 


119 


0.07 


-50,-2000 


0.237 


< 13268 


4.4 [3-7] 


1.8 [1.1-2.3] 


0.04, 70.43 


0.02, 


0.70 


120 


0.01 




0.038 


< 265 00 


1.3 [0.9-1.9] 


1.0 [0.6-1.3] 








123 


0.16 


-425,-2000 


0.585 


< 295 20 


4.7 [3-7] 


4.0 [2.6-5.3] 


7.27,161.03 


0.34, 


1.62 


125 


1.41 


-250,-2250 


5.03 


78298 


8.4 [6-14] 


19.4 [12.3-25.5] 


12.03,974.77 


0.96, 


8.64 


129 


0.37 


-425,-2250 


1.33 


87715 


4.1 [3-6] 


10.6 [5.3-14.1] 


18.98,531.97 


0.89, 


4.71 


130 


0.005 




0.017 


< 18220 


1.0 [0.7-1.5] 


0.5 [0.3-0.7] 








131 


0.24 


-175,-2250 


0.867 


62899 


3.9 [3-6] 


7.2 [4.4-9.7] 


2.20,363.81 


0.26, 


3.28 


132 


0.50 


-800,-2250 


1.79 


< 37223 


7.3 [6-11] 


7.9 [5.3-10.6] 


50.43,398.88 


1.26, 


3.54 


133 


0.12 


-100,-2250 


0.440 


<20964 


4.8 [3-7] 


3.0 [1.8-4.0] 


0.30,151.53 


0.06, 


1.39 


134 


0.58 


-700,-2250 


2.076 


75733 


5.6 [4-9] 


12.2 [7.0-15.8] 


59.63,616.12 


1.71, 


5.49 


136 


0.02 




0.088 


142440 


0.8 [0.6-1.3] 


3.5 [2-4.5] 








137 


0.10 


-800,-2250 


0.368 


< 18399 


4.7 [3-8] 


2.6 [1.6-3.3] 


16.25,128.55 


0.40, 


1.14 


138 


0.02 




0.077 


< 13800 


2.5 [1.9-3.6] 


1.0 [0.6-1.4] 








141 


0.13 


-275,-2250 


0.459 


50491 


3.2 [2-5] 


4.7 [2.6-6.2] 


3.51,235.14 


0.26, 


2.09 


142 


0.19 


-250,-2250 


0.688 


<20964 


6.0 [5-9] 


3.7 [2.5-4.8] 


2.30,186.11 


0.18, 


1.66 


143 


1.08 


-850,-2250 


3.87 


65891 


8.1 [6-12] 


15.6 [10.6-20.2] 


111.96,784.51 


2.63, 


6.97 


144 


0.09 




0.303 


<27050 


3.5 [2.5-5.5] 


2.8 [2-3.5] 








145 


0.03 


-200,-2250 


0.099 


<20964 


2.3 [2-3] 


1.4 [0.9-1.9] 


0.56, 71.31 


0.05, 


0.59 


149 


0.03 




0.094 


38090 


1.7 [1.3-2.5] 


1.8 [1.2-2.5] 








150 


0.14 


-500,-2250 


0.515 


<24816 


4.8 [3-7] 


3.5 [2.2-4.8] 


8.75,177.20 


0.35, 


1.58 


151 


1.84 


-850,-2250 


6.56 


101835 


8.4 [6-14] 


25.2 [16.7-32.6] 


180.91,1267.6 


4.27, 


11.29 


152 


0.02 


-350,-2250 


0.081 


<6417 


3.7 [3-6] 


0.7 [0.4-1.0] 


0.85, 35.29 


0.06, 


0.37 


153 


1.47 


700,2250 


5.259 


204523 


5.3 [3-8] 


32.0 [17.6-44.0] 


155.72,1608.9 


4.45, 


14.29 


156 


0.12 


-750,-2250 


0.443 


< 12408 


6.2 [5-10] 


2.3 [1.4-3.1] 


12.80,115.20 


0.34, 


1.03 


158 


0.09 


-750,-2250 


0.330 


< 15407 


4.9 [4-7] 


2.2 [1.4-2.9] 


12.32,110.84 


0.33, 


0.98 


159 


0.53 


-700,-2250 


1.90 


73980 


5.3 [4-8] 


11.6 [7-16] 








160 


0.04 


-750,-2500 


0.158 


<8130 


4.7 [3-7] 


1.1 [0.7-1.5] 


6.40, 71.11 


0.17, 


0.56 


161 


0.004 




0.015 


<1100 


3.9 [3-6] 


0.1 [0.1-0.2] 








163 


0.198 


-450,-2500 


0.706 


45779 


4.1 [3-6] 


5.5 [3.5-7.0] 


11.17,344.80 


0.49, 


2.74 
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ID 




x 10 15 


Space Velocity" 


L Ha x 10- 37 


Volume b 




MassxlO~ 3 d 


KExl0~ 51 e 


Momentum X 10 -45 e 




(erg s" 


- 1 cm- 2 ) 


(kms- 1 ) 


(erg s _1 ) 


(pc- 3 ) 


(1/V7cm- 3 ) 


(V7M ) 


(V7erg) 


(VJ dyne s) 


166 




0.271 


-250,-2500 


0.966 


62899 


4.1 [3-6] 


7.6 [4.4-10.6] 


4.71,470.98 


0.38, 3.79 


167 




2.527 


-850,-2500 


9.02 


117661 


9.2 [7-14] 


31.7 [17.6-44.0] 


227.72,1969.9 


5.36, 15.77 


169 




0.01 




0.035 


23740 


1.3 [1-2] 


0.9 [0.6-1.2] 






173 




0.089 


-750,-2500 


0.316 


<26102 


3.6 [3-6] 


2.8 [0.5-1.1] 


15.76,175.09 


0.42, 1.41 


174 




0.0003 




0.001 


<13250 


0.3 [0.2-0.45] 


0.1 [0.1-0.2] 






179 




0.089 


-400,-2500 


0.310 


<29()94 


3.4 [3-5] 


2.9 [2.2-3.5] 


4.62,180.54 


0.23, 1.43 


180 




0.027 


-275,-2500 


0.096 


36369 


1.7 [1-3] 


1.8 [1.3-2.6] 


1.39,114.95 


0.11,0.96 


181 




0.0007 




0.003 


<13250 


0.5 [0.4-0.8] 


0.2 [0.1-0.2] 






184 




0.18 


-750,-2750 


0.656 


124938 


2.4 [1.9-4] 


8.8 [5.3-12.3] 


49.25,662.09 


1.32, 4.84 


186 




0.031 


-600,-2750 


0.111 


<9843 


3.5 [3-6] 


1.1 [0.6-1.4] 


3.79, 79.52 


0.12,0.56 


188 




0.12 




0.427 


<17120 


5.3 [4-8] 


2.6 [1.5-4] 






189 




0.01 


-500,-2750 


0.037 


<24389 


1.2 [1-2] 


0.9 [0.5-1.2] 


2.19, 66.31 


0.09, 0.48 


190 




0.0001 


-250,-2750 


<0.001 


< 11548 


0.1 [0-0.2] 


0.1 [0.0-0.2] 


0.05, 6.39 


0.01,0.10 


192 




0.0008 


-750,-2750 


0.003 


38084 


0.3 [0.1-0.4] 


0.4 [0.2-0.4] 


1.97, 26.51 


0.05,0.19 


194 




0.068 


-750,-2750 


0.241 


< 10269 


5.1 [3-8] 


1.5 [1.1-2.2] 


8.37,112.56 


0.22, 0.81 


201 




0.111 


-750,-2750 


0.395 


< 30806 


3.7 [3-6] 


3.4 [1.8-5.3] 


19.21,258.25 


0.51, 1.87 


203 




0.0070 




0.025 


<22080 


1.1 [0.8-1.8] 


0.7 [0.5-1] 






205 




0.0165 


-400,-2750 


0.059 


<18833 


1.8 [1-3] 


1.1 [0.9-1.8] 


1.68, 79.47 


0.09, 0.61 


206 




0.068 


-700,-2750 


0.241 


< 10269 


5.1 [3-8] 


1.5 [1.1-2.2] 


7.29,112.50 


0.21,0.83 


207 




0.021 


-500,-2750 


0.076 


<7277 


3.4 [2-5] 


0.7 [0.4-1.3] 


1.75, 52.99 


0.07, 0.39 


208 




0.0007 


-400,-3000 


0.002 


<17112 


0.5 [0-1] 


0.2 [0.1-0.2] 


0.28, 15.85 


0.02, 0.13 


209 




0.0007 




0.003 


<12150 


0.5 [0.4-0.8] 


0.2 [0.14-0.28] 






210 




0.027 


-600,-3000 


0.097 


<7277 


3.9 [3-6] 


0.8 [0.5-1.1] 


2.83, 70.87 


0.10, 0.48 


211 




0.111 


-600,-3000 


0.395 


< 30806 


3.7 [3-6] 


3.4 [1.8-5.3] 


12.29,307.24 


0.41, 2.07 


213 




0.086 


-750,-3000 


0.307 


148475 


1.5 [1-2] 


6.6 [4.4-8.8] 


36.94,591.03 


0.99, 3.94 


214 




0.071 


-500,-3250 


0.251 


<7704 


6.0 [5-9] 


1.3 [0.9-2.6] 


3.28,138.74 


0.13,0.86 


216 




0.010 


-500,-3250 


0.036 


< 11982 


1.8 [1-3] 


0.6 [0.3-1.1] 


1.53, 64.72 


0.06, 0.40 


217 




0.114 


-500,-3250 


0.406 


<24816 


4.3 [3-6] 


3.1 [1.8-4.4] 


7.66,323.59 


0.31,2.00 


219 




0.005 




0.019 


< 10490 


1.4 [1.1-2.4] 


0.4 [0.2-0.8] 






220 




0.001 




0.004 


<9390 


0.7 [0.5-1.1] 


0.2 [0.1-0.4] 






223 




0.002 




0.008 


<7730 


1.1 [0.8-1.6] 


0.2 [0.1-0.4] 






224 




0.0001 


-400,-3500 


<0.001 


<27381 


0.1 [0-0] 


0.1 [0.0-0.2] 


0.14, 10.78 


0.01,0.08 


225 




0.013 


-450,-3500 


0.047 


<24816 


1.5 [1-2] 


1.1 [0.5-2.1] 


2.13,128.88 


0.10,0.75 


226 




0.0001 


-450,-3500 


<0.001 


<11548 


0.1 [0-0.2] 


0.1 [0.0-0.2] 


0.18, 10.65 


0.01,0.07 


227 




0.0016 


-400,-3500 


0.006 


< 15407 


0.7 [0-1] 


0.3 [0.1-0.5] 


0.42, 32.35 


0.02, 0.15 


228 




0.0001 


-400,-3500 


<0.001 


< 11548 


0.1 [0-0.2] 


0.1 [0.0-0.2] 


0.14, 10.78 


0.01,0.08 


229 




0.004 


-400,-3500 


0.015 


< 18825 


0.9 [0.7-1] 


0.5 [0.3-1.1] 


0.84, 64.03 


0.04, 0.39 


232 




0.0001 


-400,-3500 


<0.000 


<6417 


0.2 [0-0.4] 


0.1 [0.0-0.2] 


0.14, 10.78 


0.01,0.08 


233 




0.0002 




0.001 


< 18770 


0.2 [0.1-0.4] 


0.1 [0.1-0.4] 







a Space velocity from model in this paper, and paper II, respectively. Velocities were assigned from radial velocities that were mapped from the ground at 
1" resolution (see Paper II). 

b < indicates that WFPC2 did not resolve one or both of the cloud dimensions. The value tabulated was established by fitting the intensity profiles of the 
brightest filaments, or estimated by eye for the rest of the cloud population. 

Ionized gas density in units shown, mean value and [upper-lower] range shown. 

d Ionized gas mass in units shown, mean value and [upper-lower] range shown. 

e Using space velocities from this paper, and paper II, respectively. 

f This is the main filament associated with the jet. Flux listed is sum of low- and high-velocity line components, see §3.10. Velocity is the flux-averaged 
value. 
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Table 4 

Derived Parameters of the Ionized Gas Associated with the Jet 



Parameter 


Value 


Note 


[NlI]AA6548+6583/Ha 


3 ±0.2 


Typical of AGN 


Velocity FWHM 


440 km s _1 


Corrected for instrumental resolution 


Ha emission 


9.6 x 10~ 15 ergs s" 1 cm" 2 arcsec" 2 


10% of superbubble Ha flux 


Ionized mass 


(l-3)xlO 5 v7M 




KE 


(0.5- 1.5) x 10 53 V7ergs 




H recombination time 


> 10 4 yrs 





Table 5 

Derived Parameters of the Superbubble Outflow from Paper II 



Parameter Expression 
Timescale from velocities 

Adiabatic expansion timescale 
of a wind-inflated bubble 

is 

KE rate into outflow 
Wind mechanical luminosity 
Input momentum rate 
Wind+ionized entrainment mass 



Note. — L w m is the kinetic/mechanical luminosity of the wind in units of 10 43 ergs s -1 , 
L ir _\\ is the IR luminosity of the starburst in units of 10 u Lq, g nuc is the fraction of this 
from the nucleus, and (3 is the fraction of the bolometric luminosity of the nuclear region 
that is radiated by stars assuming a Salpeter initial mass function (IMF) that extends to 
100 M Q . 



^ b Rb U bbi e ,k pc /y bubble yr 

Rbubbie = 0.6(LwA3tj; ubb i efi /n Q ) - 2 kpc 
Vbubbie = 40Q(L wA3 /n tl uhble 6 ) Q - 2 km s" 1 

tbubble = 6 X l0 5 Rbubble,kpc /Vbubble,1000 VI 

Ekin ~3x \0 Al Ekin,55VbubbleM)tM/Rbubble,kpc ergs S" 
L w ,43 = 3 R bubble, kpcVbubble,10Q0 n er 8 S S 
Pkin ~3x l0 32 pkinA6Vbubble,lQQo/Rbubble,kpc dyne 
M= 10MbubblejVbubble,100o/Rbubble,kpcMQ yr" 1 
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